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Abstract 

The in-in path integral of a scalar field propagating in a fixed background is formulated in a 
suitable function space. The free kinetic operator, whose inverse gives the propagators of the 
in-in perturbation theory, becomes essentially self adjoint after imposing appropriate boundary 
conditions. An explicit spectral representation is given for the scalar in the flat space and the 
standard propagators are rederived using this representation. In this way the subtle boundary 
path integral over the field configurations at the return time is handled straightforwardly. It turns 
out that not only the values of the forward (+) and the backward (-) evolving fields but also their 
time derivatives must be matched at the return time, which is mainly overlooked in the literature. 
This formulation also determines the field configurations that are included in the path integral 
uniquely. We show that some of the recently suggested instanton-like solutions corresponding to 
the stationary phases of the cosmological in-in path integrals can be rigorously identified as limits 
of sequences in the function space. 
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I. INTRODUCTION 



The path integral formulation not only provides an equivalent alternative to the canonical 
quantization procedure but also it is an essential tool in understanding nonperturbative 
aspects of quantum theories. In addition to the vast range of applications to flat space 
quantum field theories, path integrals offer a physically transparent picture of the Hawking 
radiation [I, 2] and moreover it is possible to give a nonperturbative formulation of quantum 
gravity in terms of Euclidean path integrals [3]. 

In cosmology, however, one is not interested in transition amplitudes but rather the 
vacuum expectation values of operators. In that case, the standard formulas in perturbation 
theory and the usual expressions in the path integral formulation designed to yield the 
transition amplitudes must be appropriately modified. This revised formalism is called 
in-in (or Schwinger-Keldysh or closed time) perturbation theory or correspondingly in-in 
path integral formalism [4, 5]. An in-in path integral giving the vacuum expectation value 
of an operator is different in important ways from an in-out path integral. Namely, the 
integration variables are doubled corresponding to forward and backward propagations in 
time, the boundary conditions are completely distinct and there is a subtle integration over 
the field configurations at the return time. In this new setup, even the equivalence of the 
operator and the path integral formulations in perturbation theory, which only involves 
Gaussian integrals, turns out to be a nontrivial issue [6] (see [7-11] for early applications of 
the in-in path integrals to cosmology). 

In searching for nonperturbative quantum contributions to cosmological correlations in 
the early Universe, the in-in path integral formulation is a valuable tool. In a recent paper 
[12], we apply the stationary phase approximation to the in-in path integral corresponding 
to a scalar field in a cosmological spacetime and find out that there exist nontrivial classical 
solutions contributing to the path integral, reminiscent of the instantons in Euclidean field 
theories. In this paper, we would like to sharpen our observations by formulating the in-in 
path integral as mathematically rigorous as possible. 

One way of making path integrals relatively rigorous is to make a Wick rotation the 
Euclidean signature. It is clear that this procedure is not applicable for general spacetimes 
and obviously not suitable for cosmological backgrounds where the time evolution is highly 
nontrivial. In that case, the "convergence" of the path integral can be achieved by intro- 
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ducing a small exponentially damping factor. This is the so called the ie prescription, which 
we adopt in this work. 

An in-in path integral involves three different integration variables corresponding to the 
field propagating forward in time (labeled by +), the field propagating backwards in time 
(labeled by — ) and the (boundary) field, which equals the common value of the + and the 
— fields at the return time. The integral over the boundary field has a different measure 
and is subtle to perform. In the free theory, this integral can be handled by introducing a 
delta functional, which enforces appropriate boundary conditions for the propagators after 
the Gaussian path integral [6] (see also [12] for an alternative treatment). Even in this 
case, however, one has to make sure that the operator whose inverse gives the propagators is 
symmetric, which requires additional boundary conditions. We will show that the discretized 
path integral, which is indeed the essential approach connecting the operator and the path 
integral formulations, also implies the same set of boundary conditions. This fact is mainly 
overlooked in the literature and might have important implications. 

One of the most straightforward ways of calculating a Gaussian path integral is to in- 
troduce the "momentum space" mode functions for the field variables, which factorizes the 
functional integral into multiples of simple finite dimensional integrals. This factorization is 
actually achieved by the spectral representation of the kinetic operator. As noted above, the 
in-in kinetic operator must also be symmetric to yield the propagators of the in-in perturba- 
tion theory. This suggests to introduce a suitable function space on which the in-in kinetic 
operator acts, which can be elevated to a self adjoint operator. One may then introduce a 
suitable spectral representation to calculate the Gaussian path integral. We show that this 
procedure can be used to get the in-in propagators in the flat space. In an interacting theory, 
once the correct propagators are obtained the connection with the standard perturbation 
theory in the operator formalism can easily be established. The main benefit of this new 
formulation in terms of a function space is that it uniquely determines the field configura- 
tions that are included in the path integral in a nonperturbative way, which means that one 
may now search for other ways of approximating the full path integral around perturbatively 
inaccessible configurations as in [12]. 

The plan of the paper is as follows: In the next section we consider a real scalar field 
propagating in a classical background and discuss the boundary conditions that make the 
in-in kinetic operator to be symmetric in a suitable function space. Assuming that a self 
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adjoint extension of the operator is given, we discuss how the Gaussian path integral is 
performed using a spectral representation. We also show that the discretized version of the 
in-in path integral implies the same set of boundary conditions and thus the equivalence 
of the in-in path integral formula with the operator formalism only works when these are 
imposed. In section III, we illustrate this procedure for the flat space by calculating the 
propagators using the explicit spectral representation. In the same section, we also discuss 
the scalar propagating in a cosmological spacetime and show that the instanton-like states 
found in [12] can be identified as limits of sequences in the function space. Section IV 
contains our conclusions and comments for further directions. 



II. IN-IN PATH INTEGRAL MEASURE IN A FUNCTION SPACE 

Consider a real scalar field propagating in a general curved spacetime which has the 
canonical action 

S[<t>] = ~\J d 4 xV=g [V^V0 + m 2 2 + V{<j>)] . (1) 

Assume that the spacetime can be foliated by spacelike surfaces E t corresponding to a 
global time parameter t. It is a straightforward exercise to apply the canonical quantization 
procedure to (1). For the free field (or for the field in the interaction picture) this amounts 
to introducing the Klein-Gordon inner product 

(f,9)KG = -i J dEn»f%g*, (2) 

where is the unit normal vector to E and gE is the induced volume element. The complete 
set of mode functions /ij obey [13] 

V 2 ^ - m 2 ^ = 0, (3) 

and 

{fii,fij) K G = Sij, (h*,h*)kg = Sij, {fii,fi*) K G = 0. (4) 

Stoke's theorem can be used to show that the Klein-Gordon inner product is independent 
of the choice of E. The free field operator can be expanded as 

= J^A*iOi + K a i» ( 5 ) 

i 
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where the creation and the annihilation operators satisfy 

[a i ,a]] = 5 ij . (6) 

The ground state is defined by Oj|0>= and the Fock space of states can be build by acting 
with a\ on |0>. In cosmology (4) is equivalent to the Wronskian condition imposed on the 
mode functions. 

It is well known that in a general curved spacetime the vacuum state |0> is not uniquely 
defined. There are infinitely many set of mode functions obeying (3) and (4). Since the 
mode equation (3) is second order in time derivatives, fa is uniquely determined once fa(t ) 
and n^d^fiiito) are fixed on £ <0 at some initial time t . Specifically when £ = R 3 , fa(t ) and 
n^d^faito) can be chosen to be equal to the mode functions of the flat space, i.e /ifc(to) = 



1/a/2(A; 2 + m 2 ) and n^<9 M /ifc(to) = — «a/ (k 2 + m 2 )/2. We denote the vacuum state defined 
by suitable initial conditions at time t as \0,t >. 
Let us consider the generating functional 

Z[J + ,J-] = y/)0<O,to|0,t*>j-<<M*|O,to>j+, (7) 

where f D<f>\<f>, f»><^t»| is the identity operator constructed from the field variables at some 
return time t* and the transition amplitudes are evaluated in the presence of two independent 
external sources J + and J - , which are coupled to the field variable. Differentiating (7) with 
respect to J + and J~ having different time arguments and setting J + = J~ = give various 
vacuum expectation values of the field operators in the ground state |0, to>. By introducing 
consecutive identity operators infinitesimally distributed in the time interval (to,t*) in (7), 
the following path integral representation can be given for the generating functional 

Z[J + ,.r] = J Dtj n^ + ^e< 5 ^ + ' J+1 -^^ ] )^o[0 + (to)]^r^o)], (8) 

where the vacuum wave functional are defined as ^/o[0 ± (^o)] =<0 ± (^o)|O, to>, D(f> integral 
is performed over field configurations at the constant time hypersurface and ± integrals 
are over all field configurations starting from time t and ending at time obeying 

0+(t,) = 0-(t,) = 0. (9) 
It is also possible to rewrite the path integral as 

z[j + ,j-] = I n©0 + x^(^ + ' j+] - 5 ^^^ (io) 

J t 



where now there is no restriction imposed on the integration variables, and (10) reduces to 
(8) after integrating over (f) + (U) or (/>~(t*). 

Let us now calculate the generating functional in the free theory and let us for the moment 
forget about the vacuum wave functionals and the delta functional in (10). This leaves a 
simple Gaussian path integral, which can be expressed after integration by parts as 



Z free [J + , J'] = / JJPSexp 



to 



d 4 x \f—g 



t 



$ J L$ + $ T J 



where 





L 








j+ 


L = 


—L 


, $ = 






-j- 



and 



L = -V U V + m 2 



Writing the inverse of the operator L as 



A" 



A" 
A- 



which obeys 



LA 
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5(x,x') 
6{x,x') 



the path integral can be performed using the "standard rules" , which gives 



Z free [J + , J ] = exp 



al A x d 4 x' v 7 -? jT A J J , 



to 



(11) 



(12) 



(13) 



(14) 



(15) 



(16) 



where x and x' denote two different spacetime points. From the very definition of the 
generating functional, the propagators can be determined in terms of the free field vacuum 
expectation values as 



A ++ (x,x') = i <0,t \T(f)(x)(j)(x')\0,t >, 
A (x,x') = i <O,t o |T0(a;)0(a; / )|O,to>, 
A- + (x,x') = i <0,t \(j)(x)<j)(x')\0,to>, 
A+~(x,x') = % <O,t o |0(a; , )0(a;)|O,t o >, 



(17) 



where T and T denote time and anti-time orderings, respectively. Note that the Green 
functions obey 

A + ~(x,x') = A~ + (x',x), 

A ++ {x,x') = A ++ (x',x), (18) 
A (x,x') = A—(x',x), 

which shows that A is a symmetric kernel. 

Even though the above derivation looks trivial, there are important subtleties that should 
be addressed. The existence of the off-diagonal terms in (14) begs for an explanation since 
one would naively write a diagonal Green function corresponding to a set of decoupled path 
integrals. Fortunately, this point is clarified in the Appendix of [6], which explains that the 
form of (14) is fixed by the delta functional in (10) that couples + and 0~ integrals. Indeed, 
as shown in [6] the sole effect of the delta functional is to impose 

A ++ (Q = A~ + (Q, 

A— (t„)=A+-(*„). (19) 

On the other hand, at least in t — > — oo limit the vacuum wave functionals are supposed to 
yield a suitable ie prescription such that A ++ and A propagators are given in terms of 
the time and anti-time ordered products of the free fields as in (17) (see e.g. [14]). 

However, additional conditions are still needed to make the Gaussian path integral well 
defined and to determine the Green functions uniquely. In obtaining the second order 
kinetic operator in (11), both actions >S , [0 ± ] are integrated by parts and the surface terms 
are discarded. Although in many applications t — > — oo limit is taken and the surface 
terms vanish after imposing suitable fall of conditions at the asymptotic past infinity, there 
is no reason for the surface terms to vanish at t*. The absence of the surface terms is 
crucial to have a symmetric differential operator so that the Gaussian path integral becomes 
well defined and the corresponding Green function becomes symmetric in its arguments as 
in (18). Although in a finite dimensional Gaussian integral the quadratic integrand picks 
up the symmetric part of the matrix automatically, this is no longer true for the infinite 
dimensional case. 

Since the original path integral derivation of the generating functional requires (9) to 
hold, the surface terms at t* vanish if and only if the following boundary conditions are 
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imposed 

<j> + (U) = <j>-{U), (20) 
n»d li <f> + (t*)=n' i d li <f>-(U). (21) 

The free Gaussian path integral becomes well defined after additional boundary conditions 
are implemented at to- Although the first condition (20) is already encountered in the path 
integral derivation, the second condition (21) might appear surprising and can be thought 
to ruin the equivalence of the path integral and the operator formalisms. As we will discuss 
shortly, the opposite turns out to be true, namely (21) is enforced by the equivalence. As 
noted above, (21) is a hidden assumption that is used to obtain the propagators. 

Despite the fact that the Gaussian path integral is uniquely determined after imposing 
(20) and (21), one may still search for a more direct way of calculating the integral. Let us 
recall, for example, how the propagator of the harmonic oscillator is obtained in the in-out 
formalism. The path integral representation of the generating functional is given by 

= J Dqexp Q dt^-q 2 - % ^-q 2 + jq]j . (22) 

To perform the Gaussian integral one can Fourier transform the variables as 

1 f°° 

q(t) = -^=jjEe iEt q(E), (23) 

1 r°° 

j(t) = -= dEe iEt j(E). (24) 

V 27T J-oo 

In the transformed variables, the integral becomes 



-q{E) [E 2 - m 2 ]q(-E) + ~j{-E)q{E) 



dE . (25) 



E 

Therefore, the whole path integral factorizes into multiples of one dimensional Gaussian 
integrals for each E. While this one dimensional integral over E does not converge, it can 
be shown that in and out vacuum wave functionals, which are neglected in (22), give rise 
to ie terms guaranteeing the convergence of the ^-integral. Introducing these ie terms, the 
final result can be written as (see e.g. [15]) 

It is now a straightforward exercise to express the propagator in the coordinate basis, namely 
in the t- variable, by applying the inverse Fourier transformation to (26), which should give 
the standard result. 
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Naturally, a similar procedure is expected to be applicable to calculate the in-in path 
integral. The harmonic oscillator example shows how this can be achieved. Mathematically 
speaking, the Fourier transformation, or the set of the functions e lEt / V2tt, give an explicit 
spectral representation of the kinetic operator d 2 /dt 2 in (22). There is a continuos spectrum 1 
E G (—00, 00) and the paths included in the path integral belong to the space of square 
integrable functions in the real line: q(t) G L 2 (— 00,00). This guarantees that no surface 
terms appear after integrations by parts in (22) and the operator d 2 /dt 2 becomes symmetric, 
which has a unique self adjoint extension in L 2 (— 00, 00). 

Analogously, to evaluate the free in-in path integral (11), one can introduce the Hilbert 
space of doublets of square integrable real 2 functions, L^ 2 ' 2 \t , t*), with the following inner 
product 



<$i|$ 2 >= / d 4 x v /= <7 $1 $2, 



(27) 



R 



L = 



(28) 



where R is the region bounded by the times to and t*, i.e. R = £ x (t ,t*), $1 = (0^,0! ) 
and $1 = (0^,02"). The kinetic operator 

L 
—L 

where L = — V^V M + m 2 , is a second order differential operator acting on L^ 2 ' 2 \t ,U). To 
avoid complications related to the vacuum wave functionals, let us take to = —00. In that 
case, L becomes a symmetric operator provided (20) and (21) are imposed as boundary 
conditions. Assume now that these boundary conditions ensure L to be essentially self 
adjoint operator in L^ 2 ' 2 ^— 00, t*). In that case, one can introduce a complete set of real 
(generalized) eigenfunctions 

■u + (A, x) 
w~(A, x) 



u{\) 



(29) 



satisfying 

Lu(A) = \u{\) 
and the boundary conditions (20) and (21) 

u + (t*) = u~(t*), 
n^d^u+iQ = n^d^u-iQ. 



(30) 

(31) 
(32) 



1 Actually the spectrum is given by A = — E 2 and for each eigenvalue A there correspond to two eigenvectors. 

2 It is straightforward to construct the Hilbert space for complex functions. Since we are considering a real 

scalar field, we prefer to work with a real Hilbert space. To clarify the notation, (2, 2) stands for square 
integrable doublets of functions. 9 



The eigenfunctions u(X) span L { - 2 ' 2 \— oo, i*) and any $ e L^ 2 ' 2 \—oo,t^) obeying (20) and 
(21) can be expanded as 

$ = y d/i(A)0(A)w(A), (33) 

where dfj,(X) denotes a measure in the set of eigenvalues. The eigenfunctions can be normal- 
ized as 

<u{X)\u{X')>=5{X-X') (34) 
so that the inner product becomes 

(35) 



(36) 



<$!$>= j d(i(\)<f)(\y. 

The eigenfunctions must also obey the completeness relation 

/dfi(X) u(X,x)u T (X,x') = 5(x, x 1 ). 

Note that while <£> denotes a doublet of functions, 0(A) is a single real function. 
Using the spectral representation, the generating functional can be expressed as 

Z free [J + , J"] = fflV^exp <$|L|$> + <$|J>V 



L>0(A) exp / d/i(X) 



-_A0(A) 2 + 0(A)J(A) 



(37) 



where J(X) =<«(A)|J >. Since the spectral functions u(X) obey (31), the delta functional 
in (10) is automatically satisfied (actually the path integral is over a more restricted set of 
of functions, which also obey (21)). 

As in the case of harmonic oscillator, the Gaussian path integral factorizes into one 
dimensional ordinary integrals, whose convergence can be fulfilled by introducing a small 
imaginary piece A — > X — ie. Evaluating the integral over 0(A) one finds 

J(A) 2 



J free 



cxp 



dfx(X) 



X — ie 



(38) 



By comparing (38) to (16) and using the completeness relation (36) the propagator can be 
expressed as 

u(X, x)u T (X, x') 



A — it 



(39) 



The vacuum wave functionals entering the path integral at to = ~ oo as in (8) are supposed 
to yield the ie term in the path integral. One therefore sees that the free in-in path integral 
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can be naturally calculated using the spectral representation of L. Unless u + (\) = or 
u~(X) = for all A, there appears off-diagonal terms in (39). 

The inner product (27) should not be confused with the Klein-Gordan inner product (2). 
Similarly, the boundary conditions (20) and (21), which are imposed for the Gaussian path 
integral and depend on the arbitrary return time i*, are not directly related to the boundary 
conditions obeyed by the canonical fields in the quantization procedure. Both the inner 
product (27) and the conditions (20) and (21) are implemented to have a mathematically 
well defined Gaussian path integral. 

Till now, we only consider the free theory and discuss how one can obtain the in-in 
propagators in an alternative but mathematically more rigorous way. In the interacting 
theory, in case when perturbation theory is applicable, the set of fields that are integrated 
out should still be same so that the Gaussian integral yields the correct in-in propagators 
consistent with the operator formalism. When the theory becomes strongly coupled, even 
though the perturbation theory essentially breaks down, one can still imagine that the exact 
results must match the sum of the infinite perturbative series. In other words, one would 
expect the set of fields that are integrated out to be independent of whether the theory is 
weakly or strongly coupled (of course the value of the path integral depends on the coupling 
constants of the theory and weakly and strongly coupled limits are expected to give very 
different results.) Therefore, the exact path integral for the generating functional 

Z[J+,J-} = / JjD0 + P0-e<^ + ' J+] - s[ ^' J_1 )^ o [0+(-oo)]^*[^(-oo)], (40) 

J -oo 

must be a sum over all in-in paths that belong to L^ 2 ' 2 \— oo, £*) obeying the boundary 
conditions (20) and (21): 

(0 + ,<T) G L (2 ' 2) (-oo,t,), (f> + (Q = <f)-(U), n^d^ + {U)=n^d^-it*)- (41) 

Although the path integral measure is not defined rigorously, (41) gives a precise nonper- 
turbative criteria about the configurations included in the sum. As we will see in the next 
section, this is important in applying approximation methods in evaluating the in-in path 
integrals. 

Our discussion so far is based on the free path integral and in the last paragraph we 
made a generalization to the exact theory. The main strategy was to be consistent with the 
perturbation theory or more precisely with the operator formalism. As is well known, the 
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first boundary condition in (41), i.e. + (t*) = _ (t*), arises naturally when one derives the 
in-in path integral from the operator formalism (this can be directly seen from (7)). The 
second condition in (41), on the other hand, is merely deduced from the consistency of the 
free Gaussian path integral and we argue that it should still hold in the interacting theory. 
It is thus crucial to see how this condition arises in the full path integral, or if not how 
the equivalence of the path integral and the operator formalisms still holds. To analyze the 
issue more carefully, let us consider the discretized path integral for (40), which is indeed 
the main object connecting the operator and the path integral formalisms. In that case, 
assuming that the metric is of the form ds 2 = —dt 2 + dTi 2 and the scalar has the standard 
action (1), (40) can be written as 



t*d(pTd(t),_ e ...exp 



(42) 



where e is a discretization parameter (not to be confused with e in the ie prescription) and we 
suppress all but the time coordinate as the index. In (42) we only write the terms that depend 
on the boundary variable 0^, which are precisely the terms containing a time derivative in 
the action. Note that the terms involving the spatial gradients and the potential cancel out 
each other between ± branches at time t*. On the other hand, the source couplings turn 
out to be e-suppressed compared to (43) or one may take U large enough so that the sources 
vanish, i.e. ^(t) = for t > t*, as we assume in (42). We thus see from (42) that the 
boundary integral has the form 



/ 



# t „...exp 

which implies 



■2iK ( K ~ c ^ ) - 



(43) 



lim ^ - ^ = Inn _ ^ ~ ^-J = ^ _ ^ = Q {u) 

This is nothing but (21) written in coordinates normal to the surface S. It is remarkable 
that the path integral imposes the same conditions that are necessary to make the kinetic 
operator symmetric. 

Summarizing our results, we find that to calculate the in-in generating functional one 
must integrate over all field configurations in L^ 2 \— oo, t*) obeying (20) and (21). In per- 
turbation theory, the same conditions are required for the kinetic operator to be symmetric 
and assuming there exist a self adjoint extension, the propagators can be calculated using 
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the spectral functions. Recall that the first condition arises due to the identity operator in 
(7) and the second follows from the boundary integral as shown above. The condition (21) 
is mainly overlooked in the literature (see e.g. [16] for an exception) and it should clearly 
be considered in determining the effects of boundary terms to the cosmological correlation 
functions (see [17-21]). 



III. THE FLAT SPACE AND COSMOLOGICAL EXAMPLES 

In the previous section, we discuss how the free in-in Gaussian path integral can be 
performed using the spectral representation of the in-in kinetic operator. In this section, we 
would like to carry out an explicit computation for a scalar propagating in the flat space 
with the Lorentzian metric 

ds 2 = -dt 2 + dx 2 + dy 2 + dz 2 . (45) 

Without loss of any generality we take t* = and consider t — — 00 case. As usual, 
by applying a Fourier transformation in R 3 with the functions exp(ix.k)/(27r) 3//2 , the in-in 
kinetic operator can be transformed into 



(46) 



where u 2 = k 2 +m 2 . To our knowledge the spectral analysis of this operator in the half line 
(—oo,0) is not studied before. 

Since the spatial dependence of the functions is handled by the Fourier transformation, 
it is enough to consider a reduced Hilbert space, one for each wave number k, consisting of 
doublets of time dependent functions = (0 + (t), <j>~(t)) with the inner product 

<$ 1 |$ 2 >= / dt$l(t)Q 2 (t). (47) 

J — oo 

The boundary conditions of interest that make L symmetric are given by 

+ (O)=0-(O), (48) 
0+(O) = 0-(O). (49) 

We would like to see whether L is essentially self adjoint or not. Since L is real the upper 
and lower deficiency indices equal and it suffices to look for normalizable solutions of 

L$(A) =i$(A), (50) 
13 



which can be solved as 



<P + = ai e Zlt + a 2 e 



-Z\t 



a 3 e Z2t + a 4 e~ Z2t , 



(51) 
(52) 



where z\ = i — u 2 and z\ = —i — uj 2 . Normalizability of the solutions requires a 2 = 0,3 = 
(note that t < 0) and the boundary conditions (48) and (49) give a\ = a 4 = 0, which shows 
that there are no eigenfunctions with imaginary eigenvalues and thus L is essentially self 
adjoint. 

To determine the generalized eigenfunctions of L, one needs to solve 



Lu(A) = \u{\). 



(53) 



Writing u(X) 



(u + , u ), this implies 

d 2 u 



dt 2 
d 2 u~ 



(\-to 2 )u + , 
-(X + oo 2 )u- 



(54) 
(55) 



The solutions are given either by sin and cos functions or by the exponential. Although 
the oscillating functions are allowed, which would yield delta function normalization of the 
eigenfunctions, the exponentially growing pieces should be discarded. The solutions for 
different values of A can be found as 



A > tu 2 : u(X) 



u 2 > A > -co 2 : u{\) 



A < -u 2 : u(X) 



bi exp (tyf\ — cu 2 ) 

b 2 sin (ty/X + lo 2 ) + b 3 cos (ty/X + UJ 2 ) 

Ci sin (ty/— X + LU 2 ) + c 2 cos (W—X + UJ 2 ) 
c 3 sin {ty/ X + lu 2 ) + c 4 cos (ty/ X + tu 2 ) 



(56) 



d 1 sin (ty/—X + u 2 ) + d 2 cos (ty/—X + u 2 ) 
d 3 exp (ty/—X — to 2 ) 



The unknown coefficients must be fixed by the boundary conditions (48) and (49) together 
with the orthonormality and completeness of the eigenfunctions. As we present the impor- 
tant details in the Appendix, these requirements determine the coefficients uniquely. For 
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A > uo 2 and A < — u 2 , there is a single eigenfunction for each A: 



\>u 2 : 



A < —u : 



11(A) = 



(A + w 



2\l/4 



exp A — w 2 ) 



v 7 ^ y/^g S in (tVA + w 2 ) + cos (ty/X + oo 2 ) 



(57) 



«(A) = 



-A + o; 2 ) 1 / 4 



^/e^sin (tV-A + w 2 ) + cos (tV-A + w 2 ) 
exp A — oo 2 ) 



On the other hand, for —oo 2 < A < oo 2 there are two eigenf unctions for each A, which are 
given by 

1/4 



oo 2 > A > -w 2 : «i(A) 



oo 2 > A > -w 2 : u 2 (A) 



VA + cu 2 - V-X + uj 2 
2^A 



VA + w 2 - v/^AT 



1/2 



1/2 



sin(tv^AT^) 

„\ 1/4 

=**f) sin(tv^T^) 



2ttA 



cos (W— A + cu 2 ) 
cos (t\/A + cu 2 ) 



(58) 



As shown in the Appendix, the eigenf unctions are properly normalized with respect to the 
norm (47) 

< u(X)\u(X') >— 8(X - A') (59) 



and they are also complete 



/oo 
dX u(X, t)u T (X, t') 
■oo 



6(t - f) 
6{t - 1') 



(60) 



where in the range — oo 2 < X < oo 2 both of the eigenfunctions Ui(X) and u 2 (X) in (58) must 
be included in the integral. 

Once the complete set of eigenfunctions are obtained, the Green function can be calculated 
using (39) as 

4M = -f«*^M. (61) 

J -oc X l€ 

Similar to the harmonic oscillator example, this computation can be carried out using the 
contour integration techniques (see Appendix), which gives the in-in propagator in momen- 
tum space 



A(M') = -i 



0(t-t') iu(t-t') , 0(t'-t) -iu(t-t') 
2uj 2uj 



J_ e -iu(t-t') 
2ui 



J_ iui(t-t r ) 
2u 

0(t'-t) iu(t-t') , 0(t-t') -iu(t-t') 
2u> 2ui 



(62) 



15 



This shows that the procedure outlined in the previous section for the calculation of the 
Gaussian in-in path integrals is consistent with the operator formalism. In particular, if the 
second condition (49) would not be imposed then (62) could not be obtained. 

In principle, the above discussion can easily be generalized to cosmological spacetimes 
with the metric 

ds 2 = -dt 2 + a(t) 2 (dx 2 + dy 2 + dz 2 ). (63) 

After applying the Fourier transformation along the spatial directions, the in-in kinetic 
operator becomes 

d 2 I 3a a 
dt 2 ^ a dt 



L = 



o -(^ + fl + ^) 2 



(64) 



where to(t) 2 = k 2 /a 2 + m 2 . As before, one may define a reduced function space on which L 
acts consisting of doublets of time dependent functions $(i) = (0 + (t), 4>~{t)) with the inner 
product 

,o 

<$ 1 |$ 2 >= / dta{tf^l(t)^ 2 {t). (65) 

J — oo 

Note that a 3 term in (65) plays the role of a weight function. 

It is easy to see that if appropriate fall off conditions are imposed at the asymptotic past 
infinity, L in (64) becomes a symmetric operator provided the usual boundary conditions 
(48) and (49) are imposed at t — t* — 0. Since the mode functions are not known for general 
a(t), it is not possible to prove that (64) is essentially self adjoint by looking at the deficiency 
indices. On the other hand, for de Sitter space in the Poincare patch that has a = exp(Ht), 
the upper deficiency indices are determined by the normalizable solutions of 

^ + 3H^- + {e~ 2m k 2 + m 2 ) + = (66) 
^ + 3H^- + (e~ 2m k 2 + m 2 ) 0" = (67) 



These can be solved as 

+ = e~ 3m / 2 [ Cl HW(kr]) + c 2 Hf\k V )] , (68) 
0- = e -^t/2 [ C3H W {kr]) + c 2 Hi 2 \k V )] , (69) 

where and denote the Hankel functions of the first and the second kinds, the 
conformal time rj is given by r] = exp(—Ht)/H, \i = a/9/4 — m? / H 2 + i/H 2 and v = 
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a/9/4 — m 2 / H 2 — i/H 2 . Using the asymptotic form of the Hankel functions it is easy to 
see that the solutions are not normalizable and thus the symmetric operator (64) is indeed 
essentially self adjoint in de Sitter space. 

Till now, we have basically studied the spectral properties of the in-in kinetic operator, 
which is crucial in perturbation theory. However, as discussed at the end of the previous 
section, the formulation of the in-in path integral in a suitable function space does not rely 
on perturbation theory. We would like to illustrate how our construction can be used to 
infer nonperturbative properties, more specifically we would like to consider the stationary 
phase approximation and show how some of the instanton-like solutions constructed in [12] 
can be rigorously understood. For that, take a self interacting real scalar field with the 
potential (see Fig. 1) 

U = ^> 2 (0-0 O ) 2 (70) 

in the Poincare patch of the de Sitter space, i.e. a = exp(Ht). We focus on the zero mode 
in the path integral, which is a sum over all square integrable doublet of functions $ with 
the norm 

r° 

<$|$>= / dte mt § T (t)<&(t). (71) 

J — oo 

We would like to show that if <p c i is an arbitrary solution of the classical equations of motion 

dV 

<t> d + 3H(j) d + —- = 0, (72) 



then the in-in field 

<t>d 



(73) 



4>d 

is included in the path integral for the generating functional. This amounts to prove that 
the following integral converges for arbitrary initial data: 

rO 



/ 

J — c 



dte^tf, (74) 



It is known that (72) gives oscillations that are damped by the expansion of the Universe. 
Viewing evolution backwards in time, 4> c i tends to grow as t — > — oo, which would fight for 
the exponential factor in (74). Let us, for the moment, forget about the mass term and the 
expansion of the Universe. Then (72) becomes 

4> d + \<f> 4 d = 0, (75) 
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FIG. 1. The scalar potential V = f (f> 2 {<j) - <j) ) 2 . 

which can be solved as 

<f> cl = ASn (aJ^, -l^j, (76) 

where Sn is the Jacobi elliptic function, which periodically oscillates between — 1 and I. To 
determine the solution with the Hubble term we note that the amplitude A also determines 
the frequency of the oscillations in (76). Since we are interested in the asymptotic behavior 
of the fields as t — > — oo corresponding to large field fluctuations, i.e. C ; ~ A ^> H, the 
frequency of the oscillations becomes much larger than the expansion rate in that limit. 
Therefore, one may try an approximate solution with slowly varying amplitude A(t). Using 
such an ansatz in the equations of motion 

d + 3 #0 + = o, (77) 
one can find the following asymptotic solution 

<P cl ~ A exp(-Ht) Sn ^~ Ht ) ^ _ij , t _oo. ( 78 ) 

Including the mass term to (77) does not alter this asymptotic form since the mass term 
becomes eventually negligible compared to 4 term as t — > — oo. Eq. (78) shows that the 
integral (74) converges and thus all classical solutions are included in the in-in path integral 
in the form (73). 

In [12], the stationary phase approximation is applied to calculate the in-in path integral 
for the generating functional and it is shown that the stationary phases are given by (73) 
where C ;(— oo) equals the expectation value of the field operator in the vacuum state at 
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t = — oo. For the potential (70), there are two possible vacuum states. Assuming that one 
expands the theory around the vacuum = 0, the stationary phases must obey 

0(-oo) = 0. (79) 

As discussed in [12], the classical solutions with the initial data 

0(-oo) = 0, 0(-oo) = v (80) 

can interpolate between the vacua = and = 0o. Namely, even one is expanding around 
the vacuum at = 0, the vacuum at = O also contributes to the generating functional 
via these classical configurations. They are very similar to the instantons of Euclidean field 
theories. 

Although the construction of instanton-like solutions presented in [12] is physically viable, 
which can be attributed to the spontaneous fluctuations around the vacuum, giving initial 
data at t = — oo can be mathematically doubtful. Let us therefore analyze the problem in 
a different way. Fix an arbitrary (negative) time t and consider the classical solution with 
the initial data 

0i(t) = O, Mt)=vo, (81) 

which we identify as the first member of a series of solutions. We define the (n + l)'th 
solution by shifting the time argument of the n'th solution by 1/n: 

0„ +1 (t) = n (t-^j ■ (82) 

Since the equations of motion (72) does not depend on t explicitly, n+ i is a solution provided 
n is a solution. Since 0i is a solution by construction, all functions in the sequence must 
solve equations of motion. As proved above, all classical solutions are square integrable and 
thus belong to the Hilbert space in which the path integral is calculated. 

It is clear that as n gets larger n approaches the instanton-like configuration (80). More- 
over, 

r-0 



lim |0„ +1 - n | = lim / dte mt [0„ +1 (t) - 0„(t)f 

n— »oo n— s>oo / „ 

J — oo 



— oo 


3Ht 



lim / dt e 



1, 



4>n(t ) - 4>n{t) 



n 



0, (83) 



which shows that the sequence has an accumulation point. By completeness, the accumu- 
lation point must also be included in the Hilbert space, which shows that the in-in path 
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integral can be expanded around it. Depending on the value of the initial velocity vq, one 
has linin^oo (p n = or lim^oo <p n = <f> , which are identified in [12] as the instanton-like 
configurations for the potential (70). Namely, we see that even though for finite n, <p n does 
not obey (79) and thus the field $ = (0 n , (f> n ) is not a suitable point for stationary phase 
approximation around the vacuum at <fi = 0, the approximation gets better and better as 
n gets larger. Alternatively, in spite of the fact that = 0o does not obey the boundary 
condition (79), it is possible to find a sequence of solutions in the Hilbert space, whose 
members satisfy (79) and approach arbitrarily close to = 0o solution in the limit. 

IV. CONCLUSIONS 

In this paper, we reformulate the in-in path integral corresponding to a real scalar field 
by defining a suitable function space. This construction allows the calculation of the free 
Gaussian path integral by using the spectral representation of the in-in kinetic operator. 
Since the in-in states that are integrated out belong to a specified function space, one has 
a well defined criteria if a given configuration can be used for the expansion of the path 
integral around it. As discussed above, this is important in applying the stationary phase 
approximation to the in-in path integrals. 

An important point emphasized in this work is that there are two boundary conditions 
that must be imposed to make the in-in kinetic operator symmetric. The first condition 
(20) is obvious in the path integral derivation of the generating functional. After imposing 

(20) , the second condition (21) is also required to have a symmetric in-in kinetic operator, 
which is certainly necessary to have a well defined Gaussian path integral. We show that 

(21) naturally arises when the boundary path integral over the field configurations defined 
at the given return time is carried out. 

It would be interesting to generalize our results to vector and tensor fields propagating in 
a classical background. In that case, a function space with a positive definite inner product 
can only be introduced after removing the longitudinal components by gauge fixing. Another 
important generalization would be to include the scalar metric fluctuations, which is crucial 
in cosmological perturbation theory (see [22] for the path integral formulation of cosmological 
perturbations). There are two main technical obstacles in that case. The first one is the issue 
of the gauge freedom, which can be handled either by gauge fixing or by working with the 
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gauge invariant variables. The second one is that the exact closed form of the interactions 
are not known, i.e. one has an infinite series that can be organized in powers of the slow- 
roll parameters of inflation [23]. Therefore, first a subset of dominant interactions must be 
chosen to search for possible non-perturbative effects. Our work is in progress examining 
these issues. 

Appendix: The spectral representation and the calculation of the propagators in 
flat space 

In this Appendix we give some of the details of our computations showing that the 
functions (57) and (58) form an orthonormal set of basis vectors in L^ 2 ' 2 \— oo, 0) obeying 
the boundary conditions (48) and (49). Moreover, we prove that using these eigenfunctions 
in (61) gives the standard Green function (62). It is instructive and technically much more 
easier to study uj = case, where the operator becomes 



The eigenfunctions of are given by sia(Et), cos(Et) and exp(±Et). While sin(Et) and 
cos(Et) are allowed in the eigenfunctions of (A.l) that would yield the delta function normal- 
ization, only exponentially decaying function can show up. Since we consider the negative 
half line t G (— oo,0), one should only take the function exp(Et) with E > 0. Moreover, 
once sm(Et) or cos(Et) appear in the one entry of the eigenfunction, the other entry must 
be exp(Et) due to the minus sign in (A.l). 

Taking the above comments and the boundary conditions (48) and (49) into account, it 
is not too difficult to obtain the following eigenfunctions 



L 



dt' 2 







(A.l) 




u(E,t) 



1 sin(Et) + cos(Et) 
A exp(Et) 



v{E,t) 



1 



exp(Et) 
sin(Et) + cos(Et) 



(A.2) 



which obey 



Lu(E, t) = -E 2 u(E, t), Lv(E, t) = E 2 u(E, t), 
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(A.3) 



where as noted above we take E > 0. Using 



f cos(Et) cos(E't)dt = ^5(E + E') + ^5(E - E'), 
[ sm(Et) sm(E't)dt = ~S(E + E') + ^S(E - E'), 

J-oo 2 2 



and 



j — c 



sm(Et)dt = - — , 
E 



(A.4) 
(A.5) 

(A.6) 



which are valid for any two real numbers E and E', it is a simple exercise to prove that the 
eigenfunctions are orthonormal 



< u{E)\u{E') >=< v(E)\v(E') >= S(E - E'), 

< u(E)\v(E') >= 0, 



(A.7) 
(A.8) 



and complete 

PCX) 

/ dE [u(E, t)u T (E, t') + v(E, t)v T (E, t')] = 5(t — t'). (A.9) 
Jo 

Recall that the inner product is defined in (47). To verify (A.6), which is true in the 
distributional sense, one may introduce a small convergence factor as 

E 



sin(Et) e et dt 



e 2 + E 2 ' 



(A.10) 



If this expression is used in a convergent E'-integral, one may safely take e — > limit, which 
gives (A.6). 

When oo ^ 0, the general form of the eigenfunctions are given in (56). To determine the 
unknown coefficients in (56), one may start applying the boundary conditions (48) and (49), 
which leaves only one unknown for A > oo 2 and one for A < —oo 2 , and only two unknowns 
for oo 2 > A > — oo 2 . The normalization 



< w(A)KA') >=<5(A-A') 



(A.ll) 



fixes the eigenfunctions uniquely up to phase factors for A > oo 2 and A < —oo 2 that gives 
(57), while it imposes the following relation 



(A + oo 2 )V-\ + oo 2 + (-A + oo 2 )V\ + co 2 



a 2 + 



b 2 = -, (A.12) 

71 



V-A + oo 2 + VA + oo 2 

where a and b are two real unknowns that appear in the eigenfunctions for oo 2 > A > oo 2 as 

W\ + oo 2 sin (W-X + oo 2 ) + b cos (W-X + oo 2 ) 



u(X) 



ay/-X + w 2 sin (ty/X + oo 2 ) + 6 cos (ty/X + oo 2 ) 



(A.13) 
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FIG. 2. The contour of integration in the complex plane for f(z) in (A.14). 

In writing (A. 13) we implemented the boundary conditions. Let us remind that a and b are 
A dependent parameters. 

To proceed one should check the completeness of the eigenfunctions (60). For that it 
is useful to derive an intermediate integral identity as follows. Let f(z) be the following 
complex analytic function 



/(*) 



Z — CO z 



l\/ Z + CO 2 



2nz 



cxp (tV z - co 2 - it'Vz + co 2 ) . (A.14) 



The function has two branch points located at z — ±co 2 . By introducing two cuts as in 



Fig 2, the square roots become single valued. We define \J z ± co 2 = ■^r±exp(i6±/2), where 

r± is the distance of z to ±co 2 and 9± are measured with respect to x-axis along clockwise 

direction. It is easy to see that f(z) is analytic at z — 0. Consider now the integral of f(z) 

along the contour in Fig 2. For t, t' < no contribution comes from the large semicircle. 

Taking the imaginary part of the integral one may verify the following identity 
dX 



I 

Jul? 



2nX 



VX-cu 2 sm(t'VX + u 2 ) + VX + co 2 cos(tVA + uj 2 ) + (t <-> t') 



+ 



f 



dX 
2^A 



Vx + u 2 - \^x + 



X 



o, 



(A.15) 



cos(*V-A + to 2 ) cos(tVA + co 2 ) + sin(tV-A + co 2 ) sm(t'VX + co 2 ) 

which is true for any real co and t,t' < 0. 

Now, using (57) and (A. 13) the off-diagonal entry of the completeness relation (60) gives 

VA - co 2 sin(tVA + co 2 ) + VX + co 2 cos(tVA + co 2 )] + (t <* t') 







dX 



2nX 



+ 
+ab 



up* 

J dxi^b 2 cos^V-A + co 2 ) cos(tVA + co 2 ) + aVw 4 - A 2 sin(tV-A + co 2 ) sm{t'VX + co 2 ) 



V-A + co 2 cos(W-X + co 2 ) sin(tVA + co 2 ) + VX + co 2 sin(tv / -A + co 2 ) cos(t'V\~+ 
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Comparing this expression with (A. 15), we see that the term in the last line containing ab 
must vanish for the equation to hold. Therefore, there must exist two eigenfunctions having 
a and b terms, respectively. Using the identity (A. 15) and the condition (A. 12), which is 
demanded by normalization, it is easy to fix a and b that gives the eigenfunctions in (58). 

Since the normalization condition and the off-diagonal component of the completeness 
relation determine the eigenfunctions uniquely, the diagonal entry of the completeness re- 
lation involving 5(t — t') must be satisfied identically. The computation can be carried out 
straightforwardly although the calculation is cumbersome. Here we would like point out a 
few crucial steps. It is again useful to derive certain integral identities, which can be used 
to see cancelation of necessary terms. Let us define 

f(z) = V Z ^f exp (-i( t + t')Vz~T^) (A.16) 

and integrate f(z) in the contour shown in Fig. 3. For t + 1' < the contribution of the 
large semicircle vanishes. The differences with the previous case are that the pole at z = 
has to be avoided and the small circle gives a non-zero contribution even when it shrinks to 
zero size. Taking the imaginary part of the f(z) integral gives 

^ sm(u(t + t')) = (A.17) 

/ —V V-X + co 2 cos(yXT77(t + t')) - VX + ^ cos( v / -A + u 2 {t + if)) 
Jo 27TA L 

f°° rfAVA - CO 2 . , r- ?/ f°° + OJ 2 , rr 5/ 

~ J 2 — ^2jP\ s ™W X + " 2 (t + t'))- j 2 exp(V\-u 2 (t + t')), 

where the term in the left comes from the shrinking semicircle around z — 0. One may 
worry about the convergence of the integral for large A, but it is enough for this identity to 
hold in the distributional sense. Note that as u — > 0, (A.17) becomes (A. 6). 

Another integral identity that we use in verifying the diagonal entry of the completeness 
relation is derived as follows. Define 

2 

f(z) = exp (-i(t + t')Vz~T^) , (A.18) 

which will be integrated along the contour shown in Fig. 4. Note that there is a single branch 
point and a single cut. As before, for t + 1' < the large semicircle does not contribute. 
However, one must calculate the integral along the small semicircle since its contribution is 
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FIG. 3. The contour of integration in the complex plane for f(z) in (A. 16). 
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FIG. 4. The contour of integration in the complex plane for f(z) in (A. 18). 



finite as it shrinks down. Taking the real part of that integral gives 



Wt + t))= f 

Jo 

-f 

Jo 



dx 



UJ 



2%W\ + to 2 



cos(v / A + w 2 (t + t / )) 



dX 



00 



US' 



.cos(V-A + u 2 (t + t')). 



(A.19) 



Using the identities (A. 17) and (A.19), the diagonal entry of the completeness relation can 
be verified. 

Finally let us present some of the details of the computation that gives (62) from (61). 
This turns out to be very similar to the verification of the completeness relation (60); the 
crucial difference is the A — ie term in the denominator of (61). Let us again start from 
the off-diagonal entry of the Green function. We first observe that for an arbitrary function 
/(A) one has 

A 



>im r dAiW 



d\f(X) 



X e 
+ i 



A 2 + e 



A 2 + e 2 



/oo 
dXf(X) 
-oo 



A 2 + e 



2 

(A.20) 
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FIG. 5. The contour of integration in the complex plane for f(z) in (A. 22). 

Using this equality for A + ~ in (61), the eigenfunctions (57) and (58) give 

Im(A+-) = cos(w(t - t')). (A.21) 
To calculate the real part of A + ~, one may derive an identity similar to (A. 15). Let us define 
M = ^E?Zg^i^±g exp ( tv ^37j2 _ it'^TuA , (A.22) 

and integrate f(z) along the contour shown in Fig. 5. Note that f(z) in (A.22) is equal to 
(A. 14) divided by z — ie that generates a simple pole at z — ie. Taking the imaginary part 
of this integral gives the integrand of the real part of A + ~ multiplied by A/(A 2 + e 2 ), which 
exactly appears in the identity (A. 20). Carefully treating the pole at z = ie gives 

Re(A+") = sin(w(t - t')). (A.23) 
2cu 

Combining with (A.21) one then finds 

A+- = -— exp(ico(t -If)). (A.24) 
2cu 

In the calculation of A ++ , one encounters the following integral 

/oo J\ 
- /y~T — 2i \ ■ \ «»(VA+^(* - t')), (A.25) 
- w 2 27TV A + w^(A — ie) 

which arises 3 after using (57) and (58) in (61). By a change of variable x = vX+cu 2 ", this 
integral becomes 



I 



dx , , ... 1 f°° exr>(x(t — t')) ,. 

— cos(x(t - t')) = — / P \ 11, (A.26) 



o n(x 2 — uj 2 — ie) 2n J_ 00 x 2 — lo 2 — ie 



At first, the integral limit appears from uj 2 to oo. One adds and subtracts to get (A.25). 
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which then implies 



where the dotted term stands for other possible contributions. Let us note that if one drops 
the (A — ie) term in (A. 25), it becomes S(t — t'), which is exactly the combination giving the 
diagonal part of the completeness relation. To see that the remaining terms in (A. 27) cancel 
each other as in the completeness relation, one can modify the identities (A. 17) and (A. 19) 
by multiplying the complex functions (A. 16) and (A. 18) by l/(z — ie). By carefully treating 
the pole at z = ie and other contributions, one may verify that all these extra terms cancel 
each other and (A. 27) becomes the final answer. 

The cautious reader would recognize that (62) is equal to the complex conjugate of the 
standard in-in propagator. This is not a problem since the propagators in the position space 
will be the same after the inverse Fourier transformation. In other words, since the correla- 
tion functions are real it does not matter if one uses A or A* as the propagator. Moreover, 
once A + ~ and A ++ are determined, A~ + and A can be fixed from the symmetries of the 
Green function. 
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